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The effect of the charge ordering on the spin-Peierls (SP) state has been examined by using 
a Peierls-Hubbard model at quarter-filling with dimerization, on-site and nearest-neighbor re- 
pulsive interactions. By taking account of the presence of dimerization, a bond distortion is 
calculated variationally with the renormalization group method based on bosonization. When 
the charge ordering appears at V = V c with increasing the nearest-neighbor interaction (V), 
the distortion exhibits a maximum due to competition between the dimerization and the charge 
ordering. It is shown that the second-order phase transition occurs from the SP state with 
the bond alternation to a mixed state with an additional component of the site alternation at 
V = V c - 
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§1. Introduction 

Recently unconventional spin-Peierls (SP) states have 
been studied extensively for a system with a quarter- 
filled band, which shows a SP state in the presence of 
a dimerization and/or a charge ordering (CO). Typ- 
ical phenomena of such a state have been found in 
quasi-one-dimensional organic conductors, (TMTTF) 2 X 
(X=PF 6 ,AsF 6 ), 1,2 ) which undergo the SP transition at 
temperature being much lower than the CO temperature. 
The onset temperature of the CO state in the TMTTF 
salts 2 ) has been identified with the temperature corre- 
sponding to the anomaly of the dielectric constant, i.e., 
the temperature at which the ferroelectric transition oc- 
curs. 3, 4 ) In the several features of these SP states, the 
present paper concerns with a coexistence and/or a com- 
petition of the SP state with the CO as maintained in the 
following NMR experiment on (TMTTF) 2 AsF 6 . 5 ) With 
increasing pressures up to 0.15 GPa, the critical temper- 
ature of the SP state, T$p, increases but that of the CO, 
Tqo , decreases where the NMR line shapes indicate a co- 
existence of the SP state and the CO state. For pressures 
larger than 0.15 GPa, Tsp decreases gradually while the 
CO is absent. 

Several SP states at quarter-filling have been studied 
theoretically for a one-dimensional extended Hubbard 
model coupled with a lattice by using the variational 
calculation of the distortion. There are various den- 
sity waves. The charge density wave (CDW) is the site- 
centered wave with a maximum on the lattice site and 
the bond order wave (BOW) is the bond-centered wave 
with the maximum on the bond between two neighboring 
lattice sites. Two kinds of SP states have been obtained 
depending on the location of the phase of the 2k F wave 
and that of the 4fcp wave where k F (= 7r/4a) is the Fermi 
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wave vector with a being a lattice distance. One of them 
is a state with 2k F BOW and Ak F BOW (state (a)) and 
the other one is a state with 2k F CDW and 4k F CDW 
(state (b)) where 4fc F BOW and 4fc F CDW represent the 
dimerization and the CO, respectively. The state (a) cor- 
responds to the state with 2k F CDW2, 2k F B0W2 and 
4k F BOW of Ung et aZ., 6 ' the Di phase of Riera and 
Poilblanc, 7 ' and BCDW of Clay et al., 8 ^ while the state 
(b) corresponds to the state with 2k F CDW1, 2k F B0W1 
and Ak F CDW of Ung et al?) and 4fc F CDW-SP of Clay 
et al. % } Although the state (a) has been obtained explic- 
itly for the large on-site repulsive interaction, 7 ) it is com- 
plicated to find the state (b). The fact, that the CO state 
relevant to the state (b) appears in the case of the large 
repulsive interaction between electrons of the nearest- 
neighbor sites, has been shown by both a numerical di- 
agonalization 9 ) and a mean- field theory. 10 ) The appear- 
ance of the CO state is understood in terms of the phase 
Hamiltonian where the CO is followed by the change of 
the sign of the commensurate potential of the umklapp 
scattering of the quarter-filled band. 11,12 ) There arises a 
competition between the CO state and the SP state (i.e., 
the state (a)). 8, 13,14 - ) A phase diagram was obtained as 
the functions of the electron-phonon couplings for intr- 
asite and intersite displacements of the lattice, 8 ) where 
the state (b) is obtained only in the presence of the in.tr- 
asite displacement. By calculating self-consistently the 
state (a) in the presence of the dimerization but with- 
out the state (b), the maximum of the bond distortion 
was obtained at the onset of the CO. 13 ) The state (b) 
for the extended Peierls-Hubbard model (i.e., only with 
the intersite electron-phonon coupling) has been first dis- 
covered by Seo et al. 14 ^ who examined the region of the 
CO state. Calculating self-consistently both the dimer- 
ization and the Peierls distortion, in which the phase of 
the Peierls distortion is taken account, they obtained a 
first order transition between the state (a) and the state 
(b), corresponding to DM+SP and CO+SP respectively 
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in their notations. 14 ) Thus the competition between the 
dimerization and the CO is expected to exhibit a rich 
variety of the SP states. 

In the present paper, treating the dimerization as an 
external field, we show the detail of the origin of the max- 
imum of Tgp (the state (a)) as a function of the nearest- 
neighbor repulsive interaction based on our preliminary 
work. 13 ) Further we investigate the maximum of Tgp at 
the onset of the CO state by considering not only the 
state (b) but also the states with the arbitrary phase of 
the Peierls distortion. 14 ) The treatment of dimerization 
as the external field may be reasonable for the study of 
the SP state of the organic conductors, TMTTF-salts, 
since the dimerization does exist even at room tempera- 
tures 15 -* and Tqo occurs at much lower temperature. In 
§2, formulation is given. The Hamiltonian is expressed in 
terms of the bosonization and the renormahzation group 
(RG) equations are derived to calculate the distortion, u, 
i.e., the order parameter of the SP state. In §3, RG flows 
are calculated for both the state without CO and that 
with CO. The distortion u is calculated as a function of 
the nearest-neighbor repulsive interaction where a phase 
of the distortion is determined to obtain the optimum u. 
In §4, discussion is given. 

§2. Formulation 

We consider a quarter-filled extended Peierls-Hubbard 
model with dimerization, given by 

JV-l 



3=0 f 



X ( C ],* C 3+l,a +h.C.) 



we define phase variables 6± (x) and (f>± (x) as 
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where the notation a = +(— ) denotes spin for | (J. 
) and p+(-) : a represents the density operator for the 
right going (left going) electrons. The density opera- 
tors p+(-),<j{<l) satisfy the boson commutation relation, 

[p±A~l) P±,v'(<l')} = ±(<7£/ 27r )<W<W- In cc f ( 3 )' e + 
and (j) + represent the charge fluctuation and the spin 

fluctuation, respectively, since ir~ 1 d8 + /dx (n~ 1 d<t>+/dx) 
expresses the charge density (spin density ). 16 ) They 
satisfy the commutation relation, [9+(x), 9-(x') } — 
[<p + (x) 7 4>-{x') } — i7rsgn(x — x') . Using phase vari- 
ables, the operator for the right going (left going) elec- 
tron, p = +(— ), is written as 16 ~ 18 ) 

=exp ipkFX 



V2 



TTOS 



+ p0_+<7(0++p0_)]/2 



(4) 



where Cj i<7 = \fa ( ij}+. a {x) + ij)_A x ) ) with x — aj and 
the lattice constant a. In eq.(4), 5+j = ir(N + j + 
JV-, T )/2, S_, t = -7r(JV+, T + JV_, T )/2, H+j = ir(N+ A + 
iV_ iT ) + ir(N +d + N_ tl )/2 , H_ a = ^(iV + , T + JV_, T ) - 
k(N + ^ + iV_^)/2 and N p>a = J dx i>\ t<T i\> v ,a ■ In terms 
of these phase variables, eq. (1) is rewritten as, 11 ' 16 ~ 19 ) 



H = H + H 1 
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where 



Uj = —ucos(irj/2 — () , 



(1) 



(2) 



and 



l 3,° 



L ],<J^3-a- 



Quantities u and Q correspond to 
an amplitude and a phase of the Peierls distortion, re- 
spectively and Cj a denotes the creation operator for an 
electron with spin a at the site j. The Xd term rep- 
resents the bond dimerization while U and V are cou- 
pling constants for the repulsive interactions of on-site 
and that of nearest-neighbor site, respectively. The last 
term expresses the elastic energy induced by the distor- 
tion where K is a spring constant. We use a coupling 
constant g(= Aa/irvpK) instead of K , where vp is the 
Fermi velocity. The electron band splits into an upper 
band and a lower band due to the dimerization Xd, where 
the lower band is linearized to study the state close to 
the Fermi point. By applying the bosonization method, 



+ 
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dx 



j/i/4 cos 46^ 



2ira 2 

— 2/1/2 sin2# + + y a cos20 + 

— (4\/2au/v F ) sin(# + + Q cos(p + 
+4(au/v F ) 2 /g] , 



(7) 



where coefficients of the nonlinear terms are given by 
K p = {l/{1 + U + AV)} 1 ' 2 , 

K a = {i/(i-u)y/ 2 , 

y l/A = A\a/2a) 2 U 2 (U - W), 
yi/2^2(xd/t)U/{l + (xd/t) 2 }, 

y* = u ■ (8) 

The quantities v p and v a are velocities given by v p = 
v F (l + U + AV) 1 / 2 and v a = v F (l - U) 1/2 where v F = 
V2ta{l - (x d /t) 2 }/{l + (xd/t) 2 } 1 / 2 , U = Ua/(7rv F ), 
V = Va/(irv F ) and A = {1 - (x d /t) 2 }/{l + (x d /t) 2 }. 
The quantity K p (K a ) expresses the degree of the charge 
(spin) fluctuation. The amplitude, j/1/2, originates in 
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the umklapp scattering of half-filling due to the dimer- 
ization while y\ u denotes that of quarter-filling obtained 
by integrating the contribution from upper band. 11 -* The 
quantity y a denotes the backward scattering. The cutoff 
parameter, a, is of the order of the lattice constant and 
is taken as a — 1.9a/-7r for quarter-filling. 11 ) In terms of 
phase variables, the order parameters of 2fcp CDW, 4kp 
CDW and 4fcp BOW are written as 



p,<7 

2 

— cos(2&i?2i + 9+) cos <f>+ , 
na 



(9) 



27r 2 a 2 



cos(4fcp x + 29+) , 



(10) 



- y 2p cos 29+ + y a cos 2(f>+ 

— y ps sin 9+ cos <p+ — y pc cos 9+ cos 4>+ 

+ (<4 + 2&)/8<?l . (14) 



where 



y ps = (4:V2au/v F ) cos( , 

Upc — (4:V2au/vF) sin £ , (15) 

and the y 2p term is induced through the RG process with 
increasing the length of the scale, i.e., ae l . By apply- 
ing a transformation, a — > a(l + dl), to coefficients of 
cq. (14), 19 ~ 22 ) one obtains RG equations as (Appendix) 

d K P {1) = [2y 2 1/4 (l) + \ (y\ /2 (l) + y 2 2p (l)) 
1 



dl 



+ 



16 



:<&(0 + y P c(0) k% 1 ) > 



(16a) 



4fcF Bow = -J2 ( x - \ ) ^,1 + f ) 
p 

xV>-p,1 (x - I) ^-p.t(*+|) 

= 2^ sin M*-i) +2 *+) • < n > 

The variable x of eqs. (9) and (10), is defined at the 
lattice site while that of eq. (11) is defined on the bond, 
i.e. x = a/2, 3a/2, • • • . Equations (9) and (11) for 9+ — ► 
tt/4 show the 2k F BOW and 4fc F BOW while eqs. (9) 
and (10) for 9+ -> tt/2 represent 2fc F CDW and 4fc F 
CDW. 

Quantities u and £ in cq. (7) are determined so as 
to minimize (H) leading to two kinds of conditions 
, (dH/du) = and (dH/dQ = 0. The condition, 
(dH/du) = 0, is written as 

V2a 



u = F = ( sin(# + + () cos </>+) 



= cos sin 9+ cos ^ + ) + sin ( / cos 6> + cos ^ + ) , (12) 

where F is calculated as the function of u and £ from 
cq. (5). Instead of the condition of (dH/dQ = 0, we 
use a condition of optimizing u with respect to £, i.e., 
du/d( = 0, which is discussed later. It is noted that F 
is expressed in terms of the original Hamiltonian as 



F = 



2N ^ 



3," 



sin 



7TJ 

2 



c + 



X ( C ],<r C j+ 1 >° + h - C ) 



(13) 



The quantity u is determined to satisfy the relation, 
(y/2a/gvp)u = F. The quantity F of the r.h.s. of 
cq. (12) is evaluated by making use of the RG method, 
in which u is taken into account as the initial value. In 
order to derive the RG equations, eq. (5) is replaced by 
an effective Hamiltonian , H cS — Ho + -ffi ff , where Hf s 
is written as 



Hf = 



Vp 

2ira 2 



1 



dx 



2/1/4 cos 49+ - y 1/2 sin26>+ 



±GAi) = -y 2 Ji) I {vUl) + yUl)) 



+ 



j l Vi/2(l)= [2 - 2K p (l) + -yi/ 4 (/)Jyi/2(0 



dl 



+ ^y PS (i)y P c{i) 



y 2p (l) = 2 - 2K p (l) - -yy 4 (l) i/2 P (0 



(i&(0-i&(0) , 



d_ 

dl 



,(0 



Jl V ^ l) 



\gai) + \v2 P (l) 



(16b) 



-y 1/4 (0= [2 - SK p (l)] Vl/i (l) 

\{yl/ 2 (i)-vl P d)) , (i6c) 



(16d) 



(16e) 



±y < ,(l) = -G <r (l)y a (l)- \(yUl) + y 2 pc {l)) ,(16f) 



^j/<t(0 y ps (0+ ^ yi/2(i)y pc (i) , (i6g) 



- 7^(0 1 y P c(i) + \y\/2{i)y P s{i) , (i6h) 



where K a = 1 + G a /2 and G a = U. The initial condi- 
tions (I = 0) for these quantities are given by eqs. (8) and 
(15), and y 2p (0) = 0. In deriving eqs. (16), v p and v a are 
replaced by vp. We note that RG equations similar to 
eqs. (16) have been derived by Yonemitsu, 22 -* who treated 
the case of the phonon with finite frequency but without 
y 2p and y pc terms. Since eqs. (16) are the result of first 
order RG, the renormalized quantity increases to infin- 
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ity for the case of the strong coupling. Then when the 
renormalizcd quantity becomes large ( e.g., 2/1/2, 2/2p, \Va-\, 
[Vps + Vpc} 1 ^ 2 ^ 2 in the present calculation ), we stop 
the renormalization for the corresponding quantity. 

Here we estimate the quantity F of eq. (12) by em- 
ploying the following response function, whose limiting 
value with the long distance gives F as found in the or- 
dered state. The response function for F, which depends 
on both the imaginary time r and the space x, is defined 

by 



R(\ri-f 2 \) = (T T F(f 1 )F(f 2 ) / 

= R s cos 2 C + R sc sin 2C + R c sin 2 C , (17) 

where F(r) — cos ( sin 9 + (r) cos (f> + (r) + sin ( cos 9 + (r) 
coscj) + (r) and r = (x,vft). Quantities R S ,R C and R sc 
are given by 

R s = ( sin 9 + (r) cos <f> + (r) sin 9 + (0) cos 4> + (0) 



R c = ^ cos 8 + (r) cos cj> + (r) cos 9 + (0) cos (j> + (0) 

R sc = R cs = sin 9 + (r) cos 4> + (r) cos 9 + (0) cos <f> + (0) 

(18) 

In a way similar to eqs. (16), eqs. (18) are calculated 
where r is replaced by aexp(^). The corresponding RG 
equations are given by 23 ' 24 ) (Appendix), 



dl 



R s (l) 



dl 



RS) 



(K p (l) + y 2p (l))9(l c -l) 
+K a (l)+y a (l)]R s (l) 

+0(l c -l)yi/2(l)Rsc(l) , 

(K p (l)-y 2p (l))9(l c -l) 

+K a (l) + y„(l)]R c (l) 
+0(lc-l)yi/ 2 (l)Rsc(l) , 



(19a) 



(19b) 



K p (l)9(l c -l) + K a (l) + y a (l) R sc (l) 



+ ^0(l c -l)y 1/2 (l)(R s (l)+R c (l)) 



(19c) 

where R s {0) = R c (0) = 1/4, R sc (0) = and 9{x) = 1(0) 
for x > (x < 0) . The quantity l c denotes a scale 
at which the coupling for the charge fluctuation devel- 
ops well. Then the region with I > l c is regarded as 
the strong coupling one where the charge fluctuation is 
frozen. We take l c = min(l c i,^ c2 ) where the present 
numerical calculation for the response function gives a 
reasonable result with a choice of 2/1/2(^1) = 1 an d 
1 2/1/4(^2) I = 2 due to the first order RG. Equation (17) 
is estimated from eqs. (19), which are calculated by us- 
ing the solution of eqs. (16). Although the quantity F 
corresponds to i? 1 / 2 (|ri — r 2 \ = 00), we estimate F from 



the value of R(l m ), i-e., 

F = R 1 / 2 (l m ) 



(20) 



where l m denotes the following characteristic length. 
When 2/cr becomes of the order of unity due to the rele- 
vant 2/ps and/or y pc , the spin gap is formed and eq. (17) 
as the function of I takes a minimum at I = l m due to 
the first order RG. Since R(l) becomes invalid for I > l m , 
we replace R(co) with R(l m ). Thus, u is obtained by the 
following procedure. First, eq. (20) is used for estimat- 
ing F with various choices of u which gives the initial 
value of RG equations. Next, u is calculated by relat- 
ing F with u through eq. (12) which is equivalent to the 
Hellman-Feynman theorem. Finally, £ is determined to 
give the maximum of u. 

Equation (14) shows that the Peierls distortion is de- 
termined by y ps (l) and y pc (l) with the optimum Thus 
it is convenient to use a linear combination of y ps and 
2/p C terms. In the following section, we discuss the SP 
state by rewriting y ps and y pc terms in cq. (14) as 



2ira 2 



dx 



-y p sin (9 + +()cos(f> + 
- y pn cos(6» + + C) cos 4> + 



Vp( 1 ) 

V?/pn(0 

where y p (0) = A^au/vF and y pn (0) = 



cos ( sin £ 
— sin ( cos C 



fy P s(i)\ 
\y P c(i)J 



(21) 



(22) 



§3. Dimerization vs. CO in the SP state 

First, we note the case where the distortion is absent 
(i.e., g = and Xd =/= 0). The competition between two 
kinds of umklapp scattering given by 2/1/2 term and 2/1/4 
term is expected with increasing V, since the locking 
of (9+) of 2/1/2 term and that of 2/1/4 term in eq. (7) 
becomes different for V > U/A. There is a critical value, 
V = V c , with a fixed U where 2/1/2 becomes relevant 
(irrelevant) for V < V c (V c < V) leading to the 1/2- 
filled state of dimerization (the 1/4- filled state of the 
CO). The boundary between these two states is shown 
by the dashed curve in Fig. 1 ; 19 > 25 ) where the dotted 
curve denotes a boundary for g = Xd = 0. It turns 
out that the region for the CO state is reduced by the 
dimerization, i.e., Xd due to the competition between 2/1/4 
and 2/1/2 terms. The case of V < V c leads to the Mott- 
Hubbard state while the case of V c < V exhibits the 
CO state. The ground state of the former state is given 
by (9+) — 7r/4 due to the relevant 2/1/2^ +00) while 
that of the latter state is given by (6* + ) — tt/2 due to 
the relevant 2/1/4 (— » —00). For both states, the charge 
excitation is gapfull and the spin excitation is gapless. 

Next, we consider the case in the limit of small distor- 
tion where the perturbational treatment of the SP state 
is applicable based on the state with g = 0. In this case, 
the lattice distortion is determined so as to obtain the 
maximum value of sin(# + +£) in eq. (21) since 2/pn(0) = 0. 
Thus for V <V C and then (6+) = tt/4, one finds ( = tt/4 
leading to 2fc F BOW and 4fc F BOW (the state (a)), while, 
for V c < V and then {9 + ) = tt/2, one finds £ = lead- 
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Fig. 1. Phase diagram of an extended Peierls-Hubbard model on 
the plane of V-U with x d = 0.1 and g = 0.1 where the solid 
curve is a boundary between the spin Peierls (SP) state with CO 
(charge ordering) and that without CO. A dashed curve (dotted) 
denotes the boundary for g — > and x d = 0.1 (g = x d — * 0). 

(a) 

1 2 3 4 5 6 



(b) 

1 2 3 4 5 6 




(c) 



1 2 3 4 5 6 




Fig. 2. Spatial variation of the SP state with £ = 7r/4 for V < V c 
(the state (a)), that with £ = for V c < V (the state (b)) 
and that with < f < n/4 for V c < V (the state (c)), respec- 
tively. Dotted line denotes the 4k f charge density wave and the 
vertical arrows (horizontal arrows) denote spin ( the lattice dis- 
placement). The symbol, X ,is the location for the center of the 
singlet state. 



ing to 2k F CDW and Ak F CDW (the state (b)). These 
two states are shown in Figs. 2(a) and 2(b), respectively 
where the arrow denotes a spin and the cross denotes 
the location for the center of the singlet state. However 
the presence of the dimerization suggests a possibility of 
a state given by Fig. 2(c) (the state (c)), i.e., a mixed 
state with < ( < it/ 4. In this section, we examine not 
only the states of Figs. 2(a) and 2(b) but also the state of 
Fig. 2(c) by calculating the optimum £, which leads to a 




/ 

Fig. 3. The I dependence of K p (l),K a {l),y 1 / i {l),y 1 / 2 {l),y2p{V), 
y a (l),y pB (l) and y pc (l) for V = 2 with U = 5,x d = 0.1 and 
g = 0.1 where ( = vr/4 (u = 0.0029) (a) and C = (u = 0.0018) 
(b). Both cases lead to the bond order. 



maximum of u and then a minimum of the total energy. 

Now we calculate numerically the distortion, u, by tak- 
ing t = 1 and a = 1. Wc choose parameters as U = 5, 
Xd = 0.1 and g = 0.1 (i.e., K ~ 9.05) and the results for 
other parameters are discussed later. Figures 3(a) and 
3(b), which are obtained from eqs. (16), show the RG 
flow (i.e., I dependence ) for the case of small V{< V c ). 
The magnitude of u is taken from the solution of eq. (12) 
where u appears only in the initial value of RG equa- 
tions. The most relevant quantity is yi/ 2 and then 
2/1/4(0 changes the sign from a negative value to a posi- 
tive value with increasing /. Figure 3(a) shows £ = 7r/4 
and (0+) = n/A due to y ps {l) = y pc (l) while Fig. 3(b) 
shows C = and (9 + ) = 7r/4. Since (6 + ) = tt/A regard- 
less the initial condition of (, both Figs. 3(a) and 3(b) 
show Ak F BOW and 2fc F BOW. By substituting y ps (l) 
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and y pc (l) of Figs. 3(a) and 3(b) into eq. (22), it turns 
out that Fig. 3(a) shows the largest y p (l) indicating the 
ground state with £ = 7r/4. The relevant 2/1/2 results in 
the charge gap determined by the dimerization, which 
leads to an effectively half-filled band. A spin gap is ob- 
tained from the relevant y a , which comes from y ps and 
y pc as seen from eq. (16f). The spin gap is much smaller 
than the charge gap since 2/1/2(0 increases rapidly com- 
pared with y p (l). 

The RG flows for the case of large V(> V c ) are shown 
in Figs. 4(a), 4(b) and 4(c) where u is chosen from the 
solution of eq. (12). In all these figures, 2/1/2(0 becomes 
irrelevant due to the relevant 2/i/4(0(< 0) indicating 4fcp 
CDW with (8+) = it/2. Figure 4(a) with C = tt/4 shows 
y ps (l) = y pc (l) which implies y p (l) > and y pn (l) = 0. 
Figure 4(b) with ( = also shows the relevant y p and the 
irrelevant y pn . From the comparison of y p of Fig. 4(a) 
and that of Fig. 4(b), it is found that Fig. 4(b) gives 
the distortion larger than Fig. 4(a). However the larger 
distortion is expected from the state shown by Fig. 2(c) 
with < ( < 7r/4, since the dimerization breaks the 
symmetry around the distortion given by the cross of 
Fig. 2(b). Actually an example of such a state is shown 
in Fig. 4(c) where y p (l) increases rapidly compared with 
that of Fig. 4(b). For Fig. 4(c), y pn (l) of eq. (22) de- 
creases from zero to a negative value, while y pn (l) — 
for Figs. 3(a) and 4(b). Thus the state for V > V c is 
given by Fig. 4(c) with ( > due to the energy gain 
from both y p term and y pn term. 

The distortion u is determined by eq. (12), with F es- 
timated from the response function R(l) where eq. (17) 
is calculated from cqs. (18) and (19). In Fig. 5, the re- 
sponse function R(l) for £ = 7r/4 is shown for the case of 
small V(< V c ) (solid curve) and large V(> V c ) (dotted 
curve), respectively, where V c — 3.2. When the response 
function enters into the region of forming the charge gap, 
we stop the RG treatment for the charge fluctuation and 
take into account only the spin fluctuation by keeping 
the relevant coupling of eq. (14) as a constant for I > l c . 
The arrow indicates a location for I — l c at which the 
charge fluctuation becomes frozen due to a formation 
of the charge gap. Finally, we stop calculating the RG 
equations of eqs. (19) when it takes a minimum at I = l m 
(end point) corresponding to the SP state. It is expected 
that the difference between R(l m ) and the correct value 
at the long distance is small since the formation of the 
spin gap occurs at I ~ l m due to the SP state. Figure 6 
shows u dependence of F(= R}/ 2 (l m )), which is obtained 
from Fig. 5 with fixed C = 7r /4 and some choices of V. 
The quantity F increases monotonically with increasing 
u where the power law as a function of u is expected in 
the presence of the interaction as seen also for the half- 
filled case. 26 ) It turns out that F takes a maximum as a 
function of V. Such a maximum originates in a fact that 
the SP state with (8+) = 7r/4 is replaced by the CO with 
(8+) = it/2 in the case of V > V c . 

Figure 7 is the main result of the present calculation 
where the solution u for eq. (12) is obtained from the 
intersection in Fig. 6. The solid curve corresponds to u 
for ( = tt/A (Fig. 2(a)) while the dotted curve denotes 
u for ( = ( Fig. 2(b)). The latter case of C = is 




Fig. 4. The I dependence of K p (l), K a (l), 3/1/4 (Z), 3/i/2(0> 2/2p(i)> 
2/ct(Z), Vps(l) an d Vpc(l) f° r V = 3.7 with U = 5,Xd = 0.1 and 
g = 0.1 where ( = tt/4 (it = 0.0012)(a), C = (it = 0.0017) (b) 
and ( = 0.6 X tt/4 (u = 0.0020) (c), respectively. In the first 
panel (a), V2p(l) = 0. All these cases lead to the CO due to 
relevant 7/1 iA — > — noV 
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as 




0.004 



u 



0.002 




Fig. 5. The / dependence of R(l) for £ = it/A with V=2 (solid 
curve), 3.7 (dotted curve) where U = 5, = 0.1 and g = 0.1. 
The arrow denotes l(= l c ), at which the charge fluctuation is 
frozen. The end point of R(l) denotes a minimum of R(l) at 



Fig. 7. The V dependence of u for f = 7r/4 (solid curve), £ = 
(dotted curve) and the optimum value of £ (dashed curve) where 
U = 5, a;^ = 0.1 and </ = 0.1. 



0.03 



0.02 



0.01 




0.005 



u 



Fig. 6. The quantity F of eq. (12), as the function of u for C, = 
tt/4 with U = 5,x d = 0.1 and V = 1,2,3 and 4. The dash- 
dotted ( dash-two dotted) line denotes ( \f2aj gvp)u of eq. (12) 
for g = 0.1 (0.2), where the solution of eq. (12) is obtained from 
the intersection. 



also calculated in a way similar to Fig. 6. Both curves 
exhibit a maximum at V = V c shown by the arrow 
where the CO state appears for V > V c . The maxi- 
mum at V = V c originates in the competition between 
the dimerization and the CO, which is comprehended as 
follows. The numerical results together with eqs. (16c) 
and (16d) show that 2/1/2 becomes irrelevant and 1/1/4 
becomes relevant for V c < V. Since the length I for the 
minimum of R(l) corresponding to u is nearly equal to I 
with \y a (l)\ — 2 , the distortion is essentially determined 



by [vUl) 



2/p C (0] m e( l- (16f). Equations (16g) and 
(16h) show that y ps (l) and y pc {l) are determined by the 
umklapp scatterings, 2/1/2(0 and y2p{l), where y^ p is in- 
duced by 2/1/4 for V c < V (eq. (16c)). Thus the sudden 
decrease of u for V > V c originates in the fact that the 
increase of y 1/2(1) for V < V c is much larger than that 
of ?/2p for V c < V, close to V = V c . The increase of 
u with V just below V c comes from the decrease of K p , 
which results in the suppression of the charge fluctuation. 
From the comparison of the solid curve with the dotted 
curve, the maximum of u for V < V c is given by the SP 
state with £ = 7r/4 (solid curve) shown in Fig. 2(a), 13 ) 
while u with ( = ( dotted curve) becomes larger than 
that of the solid curve for large V (> V c ). The first 
order transition might be expected between these two 
states if the dimerization were absent. 14 ) However, the 
mixed state shown in Fig. 2(c) is rather promising for 
V > V c since the dimcrization does exist even for the SP 
state with the CO. The dashed curve denotes u of such a 
mixed state, which is obtained by choosing the optimum 
C leading to the maximum of u. The difference between 
the dashed curve and the dotted curve becomes small for 
V^tA. With increasing V for V > V c , the distortion u 
increases again since the state of Figs. 4(b) and (c) are 
compatible with the CO, and K p keeps decreasing with 
increasing V. The minimum of u for V < V c comes from 
the change of the sign of 2/1/4(0) as the function of V. 
Thus we obtain the transition from the state of Fig. 2(a) 
to that of Fig. 2(c) at V = V c . 

In Fig. 8, the V dependence of £ for the mixed state 
(Fig. 2(c)) is shown where Xd = 0.1(solid curve) and Xd — 
0.2(dotted curve). The quantity £ begins to decrease from 
7r/4 at a critical value and reduces to zero monotonically. 
Within the numerical accuracy, the onset of the mixed 
state is given by V = V c indicating the second-order 
transition. With increasing Xd, V c increases since the 
SP state with £ = 7r/4 is enhanced by the dimerization 



8 



Muneo SUGIURA and Yoshikazu Suzumura 




°0 0.5 1 

°0 2 4 

V 

Fig. 8. The V dependence of the optimum £ for Xd = 0.1 (solid 
curve) and x^ = 0.2 (dotted curve) where U = 5 and g = 0.1. 
The inset denotes the corresponding f dependence of u for V = 
2(solid curve) and V = 3.7(dotted curve) for x^ = 0.1 where the 
arrow denotes the location for the optimum f . 




Fig. 9. The V dependence of the charge gap A p and the spin 
gap A CT for f = 7r/4 (solid curve), ( = (dotted curve) and the 
optimum ( (dashed curve) where U = 5,Xd = 0.1 and g = 0.1. 
The difference among three curves of A p is invisible. The spin 
gap, which is multiplied by 10, is similar to u in Fig. 7. 



(oc 2/1/2)- The inset shows u as a function of £ where 
the location of the optimum ( is shown by the arrow. 
For V < V c , u as a function of £ increases monotonically 
while u for V c < V takes a maximum in the interval 
range of < ( < 7r/4. 

The energy, vpa~ l exp(— l m ), corresponding to the 
SP state is nearly equal to that of the freezing of the 
spin fluctuation while the charge fluctuation is frozen 
at higher energy as seen from Fig. 5 (the location in- 
dicated by the arrow). Figure 9 shows the charge 



gap, A p , and the spin gap, A CT , which correspond to 
Fig. 7. These gaps are calculated from A p = ut c exp[— l p ] 
and A (j = w c exp[— l a ] where l p = min(Z p i,Z P 2) with 
2/1/2 (W = 2, |yi/4^p2)| = 3 and \y a {la)\ = 2. The cutoff 
is chosen as oj c = 2.2, by noting that lu c = 5.3 for the 
half-filled case 27 ) and the ratio of the bandwith of half- 
filling and that of quarter-filling is 0.41. 28 ) The length 
l p is essentially the same but is slightly larger than l c in 
Fig. 5, which is chosen for the convenience of the numer- 
ical calculation of R(l). The charge gap is much larger 
than the spin gap, which is nearly equal to u. The dip 
of the charge gap and the cusp of the spin gap, which 
are found at V — V c , originates in the fact that the fixed 
point of 2/i/4(0 changes from +00 to —00 at V = V c (or 
cquivalently from +00 to for 2/1/2 (0)- 

Finally, we discuss the relation between the method of 
obtaining the optimum £ and that of (dH/dQ = shown 
just above eq. (12). The estimation of the latter one 
for the state of Fig. 2(c) is very complicated within the 
present scheme of RG method. However, the condition 
given by (dH/d() = (i.e., (cos(6»+ + C)cos0+) = 0) 
seems to be compatible with the former result with < 
( < 7r/4 if we note the following property of the quantum 
fluctuation of 9 + for V > V c . The fluctuation with low 
energy (i.e., large I) exists around (6 + ) — n/2 due to the 
fixed point of the CO state while the fluctuation with 
high energy (i.e., small I) exists around (9 + ) = 7t/4 due 
to the dimerization as seen from y 1/2(1) ^ for small 
I in Fig. 4(c). Further, it should be noted that such 
a competiton of the dimerization and the CO becomes 
noticeable for the state with V just above V c . 

§4. Discussion 

Using an extended Peierls-Hubbard model where the 
presence of the dimerization is assumed as a model of 
organic conductors, we have examined the SP state and 
obtained the following results. With increasing nearest- 
neighbor interaction, V, the transition from the SP state 
of Fig. 2(a) into that of Fig. 2(c) occurs at V — V c cor- 
responding to the appearance of the CO where the tran- 
sition is of second-order due to the dimerization. The 
competition of the dimerization with the CO results in 
the maximum of u at V — V c and the rapid decrease of 
u for V just above V c . 

We have also examined the SP state with other pa- 
rameters. The calculation of u similar to Fig. 7 for 
U =4,5 and 6 shows that u as the function of U de- 
creases monotonically in the region around the maximum 
(i.e., V = V c ) as found for the half-filled case. 26 ) The 
dimerization Xd increases V c as seen from Fig. 1 where 
V c for Xd = 0.1 (solid curve) is much larger than that for 
x d -» (dotted curve). 9 ' n ) The SP state without the CO 
is also enhanced in the presence of g since the SP state 
of Fig. 2(a) is compatible with the Mott- Hubbard state. 
The magnitude u of Fig. 7 obtained by g = 0.1 is rea- 
sonable for the organic conductor in which t <~ 2000i4T 15 ' 
and u(~ T SP ) - 10 K leading to u/t - 0.005. ^ The 
small effect of g on V c is found from the comparison of 
the solid curve with the dashed curve in Fig. 1 and from 
the fact that the enhancement of V c is about 0.1 even for 
g = 0.5. However, u is strongly enhanced by g. The g 
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Fig. 10. The g dependence of u for V = 0, 3 and 4 with U = 5 
and xa = 0.1. The inset denotes the optimum f corresponding 
to the main figure. 



dependence of u is shown in Fig. 10 with some choices 
of V. With increasing g, u increases with a power law 
due to the presence of U while u oc exp(— 1/g) for the 
conventional Peierls state (i.e., U = V = 0). The V 
dependence of u is weaker than the g dependence of u. 
Thus the critical value of V c in Fig. 7 does not depend 
much on g indicating also the small effect of g on the 
boundary in Fig. 1. 

We comment on our treatment of the first order RG. 
We obtained the cusp of u in Fig. 7 and the dip of A p in 
Fig. 9. Such an anomaly originates in the fact that the 
2/i/2 term competes with the yi/4 term. This competi- 
tion has a common feature with the Ising transition, 29 ' 30 ) 
which leads to the vanishing of A p at V = V c . Thus it is 
expected that u in Fig. 7 remains finite but the tangent 
becomes infinity. 

Finally, based on the result of Fig. 7 and Fig. 9, we 
comment on the SP state of organic conductors with the 
dimerization by noting that u oc Tsp, and pressure de- 
creases V/t due to the increase of t. The experiment 
on (TMTTF) 2 AsF 6 5 ) is interpreted as follows. The SP 
state under low pressure ( < 0.15 GPa), which shows the 
increase of Tsp and the decrease of Tqo with increasing 
pressure, corresponds to the mixed state of Fig. 2(c) ob- 
tained for V just above V c . The SP state under high 
pressures ( > 0.15 GPa ) showing the decrease of Tsp as 
the function of pressure, corresponds to the SP state of 
Fig. 2(a). Further we note the SP state of (TMTTF) 2 X 
with X = PF 6 and AsF 6 have T co ^ 62K and 103K and 
Tsp ~ 18K and 11K, respectively at ambient pressure. 5 ) 
These two SP states may be described by the mixed state 
of Fig. 2(c) since V/t of (TMTTF) 2 AsF 6 is larger than 
(TMTTF) 2 PF 6 due to the effective pressure. 
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Appendix: Renormalization Group Equations 

First, we derive the renormalization group (RG) equa- 
tions for H cS = H n + Hf s (eq.(14)), using the response 
function 21 ) given by 

R{ri - r 2 ) = /T T e w +^e~ ie + ( ^ 



H otf 

T Te P+(.i\) e -i0+fr) Sl \ (A4) 
(d/)o \ ' 

T T is time ordering operator, f = (x,vpT) and ( ) de- 
notes the average over Ho . In the absence of H° , R(fi — 
r 2 ) = exp{-K p U(ri-f2)} with U(r) = ln(\r\/a) . Si = 
T T exp[- / d 2 f-Hf] with Hf = (v F /2Tia 2 ) J dxHf and 
d 2 f — dxvFdr /2na 2 . 

We expand the nonlinear terms of treating as per- 
turbation and rewrite the response function as R(f\ — 
r 2 ) = cxp{— K c ^U{r\ — r*2)}- K° s is written by 
2/1/4)2/1/2) Assuming the scale invariance of response 
function, one obtain the RG equations. 

By expanding Hf 3 up to third order, the response 
function is written as 



i + Wrf 

e=± J 



i?(ri-r 2 )=cxp(-./p(l,2)) 

x{iy 2 /4 exp(-16/ p (3,4)){ cxp(4e[/ p (l,3) 
-/ p (l,4)-/ p (2,3) + / p (2,4)])-l} 

+ \(yi/2 + vl P ) exp (-4/ P (3, 4)) { exp (2e [/ p (l, 3) 
-/ p (1,4)-/ p (2,3) + ./p(2,4)])-1} 

+ ^(y 2 P s + Vic) cxp { -/p(3, 4) - / ff (3, 4) } 
x{ cxp(e[/ p (l,3)-/ p (l,4) 
-/ p (2,3) + / p (2,4)])-l}} 

+ J d2 hd 2 r A d 2 ? b | ^(2/1/42/1/2 _ Vi/iVlp) 

x exp ( -8.f p (3, 4) - 8.f p (3, 5) + 4/ p (4, 5) ) 
x{ exp(2e[2/ p (l,3)-/ p (l,4)-/ p (l,5) 
-2/ p (2,3) + / p (2,4) + / p (2,5)]) -1} 

- -^{V2 P y 2 ps - V2pV 2 pc - 2y 1/2 y ps y pc ) 

x exp ( -2.f p (3, 4) - 2.f p (3, 5) + / p (4, 5) - / ff (4, 5) ) 
x{ exp(e[2/ p (l,3)-/ p (l,4)-/ p (l,5) 
-2/ p (2,3) + / p (2,4) + / p (2,5)])-l} 

- ^(VvVps + y*ypc) 

x exp ( -/ p (4, 5) - 2/ CT (3, 4) - 2/ ff (3, 5) + / CT (4, 5) ) 
x{ exp(e[/ p (l,4)-/ p (l,5) 
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-/„(2,4)-/„(2,5)])-l}} 



(A-7) 



(A-2) 



where f p („)(i, j) — K p ^ a )U(fi — fj). Noting the renor- 
malization form, one obtain 

eff 1 , 2 2 \ f°° dr f r \ i+4AT P 

vi/4 - 2/1/4 + 4 (2/1/2 - y 2p ) J a -{-) 



off 



2/1/2 = 2/1/2 + p 2/1/4 2/1/2 



+ ^ VpsUpc 



dr / r\ 
a \aJ 

2/ 2p =2/2„- ^ 2/1/4 2/2^ -(-j 



dr /r\!- 8K p 
a Va/ 

i+k p -k ct 



-^(Vps-Vpc) I —(- 



dr ( r\ i+K P -K a 
a \a, 



yf s =y P s + \ (2/1/2 2/pc - 2/2 P 2/ps) / — (-) 
1 Z" 00 /r\i-2^ 



1-2K„ 



' 2/tT 2/ps 



dr f r\ 
a Va/ 



ape 



Vpc + 2 (2/1/2 2/p S + y2 P y pc ) 

1 f°° dr /r\ 



(A-3) 



By reexponentiating the renormalized form of eq. (A-2), 
one obtain 



R(r 1 -r 2 ) = e X p(-KfU(r 1 -f 2 )) , (A-4) 



where 



Kf = K n 



' 2/1/4 



3-16 K„ 



00 dr /r \3-4K p 
a Va/ 



- 2 (Wl/2 + V2 P ) 

- V6 (yp S + y pc ) J -(-J 



3-K„-^ 



if 2 

P 



(A5) 



From the condition of scale invariancc for a — > a(l + d/), 
we obtain the RG equations of K p , K a , y\u, 2/1/2, 2/2 P > 2/ps 
and j/pc in eqs. (16). RG equations for K a and y a are 
obtained in a similar way from the response function, 
R(ri - r 2 ) = (TT-e^+^^e-^^ 2 ) > J?eff . 

Next, we derive the RG equation of response func- 
tion 23 ) R s written as 

Rs(n - r 2 ) 

= (T r sin^ + (ri)cos^ + (fi)sin^ + (r 2 )cos(/> + (r2)) J jeff • 

(A-6) 

We calculate R s up to the second order in Hf 3 . In a way 
similar to RG equations, one obtain 

R s (n - r 2 ) = l - exp (-f p (l, 2)-/„(l, 2)) x £ a (fi - r 2 ) , 



where 

4(ri-r 2 ) = l-i y d 2 f 3 
x { j/ 2p exp ( 2 [/ p (l, 2) - / p (l, 3) - jp(2, 3) ] ) 
+ VrT exp ( 2 [f a (l, 2) - /„(!, 3) - /„(2, 3) ] ) } 



3-16if„ 



3-4K 



. 2 f°° dr /r 

4 ( ^ 2+ ^r?(D 

+^(y 2 ps + y 2 PC ) I -(-) }* 

f 1 2 f 00 dr / r \ 
x C/(fi - f 2 ) . 



3-4JC, 



. 3-K„-K a 



(A-8) 



The scale transformation, a — > a' = a(l + dZ), leads to 
(r = |ri — r 2 | = ae l ) 

R s {r,a) = I(dl)R s (r,a r ) , (A-9) 



where 

I(dl) = exp 



-y 2p dl ~ y a dl 



+ { [ 2?/i/ 4 + g (2/1/2 + 2/2 2 p) + ^(2/p 2 s + 2/p C ) ) A 2 



r ( ^2/ 2 + ^(2/ 2 s +2/ p 2 c )) A 2 K< r ' - 



From eqs. (A-7), (A-9) and (16), we obtain 

Rs(r) = 
1 



(A-10) 



■ exp 



ln(r/a) 

dl {K p {l) + K a (l)+y 2p {l)+ yry {l)} 



(All) 



which leads to 



-R s (l) =-[K p (l) + K a (l)+y 2p (l)+yAl)]Rs(l) ■ 

(A-12) 

We note that the additional term coupled to R sc can be 
calculated in terms of the operator product expansion. 24 ) 
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